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Abstract 



We calculate the non-forward quark matrix elements of operators with two co- 
variant derivatives needed for the renormalisation of the second moment of gen- 
eralised parton distributions in one-loop lattice perturbation theory using Wilson 
fermions. For some representations of the hypercubic group commonly used in 
simulations we determine the sets of all possible mixing operators. For those repre- 
sentations the one-loop mixing matrices of renormalisation factors are found. Due 
to non-vanishing contributions of operators with external ordinary derivatives the 
number of contributing operators increases compared to forward matrix elements. 
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1 Introduction 



In recent years generalised parton distributions (GPDs) £Q became a focus of both ex- 
perimental and theoretical studies in hadron physics. For an extensive up-to-date review 
including a comprehensive list of references see |2j. GPDs provide a universal (with 
the meaning used in factorisation proofs), unifying parametrisation for a large class of 
hadronic correlators, including e.g. form factors and the ordinary parton distribution 
functions. Thus they provide a solid formal basis to connect information from various 
inclusive, semi- inclusive and exclusive reactions in an efficient, unambiguous manner. 
Furthermore they give access to physical quantities which cannot be directly determined 
in experiments, like e.g. the orbital angular momentum of quarks and gluons in a nucleon 
(for a chosen specific scheme) and the spatial distribution of the energy or spin density 
of a fast moving hadron in the transverse plane. This enormous potential motivates the 
ongoing dedicated investigation of exclusive reactions at DESY, CERN, JLab and other 
accelerator centers [SI El- As GPDs are well-defined QCD objects it was possible to derive 
many fundamental theoretical results, e.g. the form of their NLO-Q 2 -evolution equations 
as well as the NLO coefficient functions for Deeply Virtual Compton Scattering. 

However, the direct experimental access to GPDs beyond the limiting cases of distri- 
bution functions and simple form factors is limited. So one can at most hope to compare 
experimental data with suitable parametrisations or models. But even under the most 
optimistic assumptions one would certainly need of the order of 20 parameters per flavour 
for a reliable fit of all GPDs. (Basically, because GPDs contain so much physics they 
cannot be expected to be trivial functions.) Another practical problem is that exclusive 
cross sections typically fall so rapidly with Q 2 that only moderate Q 2 values can be stud- 
ied. For these, however, higher-twist contributions can be sizeable, which complicates the 
situation even further. 

Therefore, although there is in principle an enormous number of reaction channels 
which provide information on GPDs, in practice it is indispensable to obtain comple- 
mentary information, e.g. from lattice QCD calculations. On the lattice we can compute 
matrix elements of local composite operators, and moments of GPDs can be related to 
such matrix elements taken between states of different nucleon momenta and spins. 
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More precisely we can write [3] for example 

i=0 

x [V] 

i=0 

+ C n (t)Mod(n + 1, 2)^u{p')u(p) A (m • • • A Mn) , 

where all indices are symmetrised and trace terms are subtracted as indicated by (■•■)• 
The leading twist-two operators used in ((TJ) are 

(. \ n-l 
\ ) $70* D» 2 ■ ■ ■ 5^) (2) 

with the symmetric covariant derivative 

3=3 - B . (3) 

Analogous equations exist for 

• \ n-l 

2 



and for the tower of operators involving a^ v (related to the generalised transversity) jE] . 
In ([T]) we have the nucleon 4- momentum transfer A = p' — p and its invariant t = A 2 , 
p = (p' + p)/2 denotes the average nucleon momentum and M its mass. 

The generalised form factors A nt2 i{t), B n ^i{t) and C n (t) are related to the moments 
of the GPDs by 1 

/l Y2~\ 
dxx^Hfa&t) = J2A nj2i (t)(-2Z) 2i + Mod(n+l,2)C n (t)(-2£) n , 
1 i=0 

dxx^Efa&t) = Yl B n , 2l (t)(-20 2l -Mod(n + l,2)C n (t)(-20 n - (5) 

i=0 



1 As an example we give here the off-forward parton distributions H (x, £, t) and E(x, £, i), one of which 
reduces to the ordinary quark distribution in the limit £ — > and t — > 0: -ff (a;, 0, 0) = 
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The variable £ = — n • A/2 is defined with the help of a light-like vector n M which obeys 
p ■ n = 1. 

First results for moments of GPDs obtained on the lattice were published recently j?HH] 
(see also and soon results from improved calculations should become available. Not 
surprisingly, a number of theoretical problems still has to be settled. One urgent task 
is to obtain the missing renormalisation factors and to completely analyse the operator 
mixing. The renormalisation of the operators which are related to (generalised) parton 
distributions has been discussed extensively, both in the continuum and on the lattice. 
Up to now, however, almost exclusively the case of forward matrix elements has been 
considered. When non-forward matrix elements are studied, new features arise, which 
make a reconsideration of the renormalisation problem necessary. In particular, the mixing 
with "external ordinary derivatives", i.e. with operators of the form d^d u • • • \M) ■ • • ip\, 
needs to be investigated JJJ. They do not contribute in forward matrix elements, but 
have to be taken into account when calculating e.g. the generalised form factors A n>2 i(t) 
for % > 0. 

On the lattice the mixing patterns are usually more complicated than in the continuum, 
because covariance under the hypercubic group H(4) imposes less stringent restrictions 
than 0(4) covariance. The necessity to consider also operators with external ordinary 
derivatives enlarges the set of contributing operators even further. These complications 
do not yet arise for n — 1 and n = 2. Hence for these moments the renormalisation factors 
can be taken over from the forward case. 

In this paper we investigate the renormalisation problem for n = 3 within the frame- 
work of one-loop lattice perturbation theory. Some first results have been presented 
recently [TO] - We find that the numerical values of the renormalisation factors are not 
very different from one. We also find that the pattern of operator mixing is far more 
involved than e.g. for moments of the ordinary parton distributions. While the numerical 
results from one-loop lattice perturbation theory have limited precision, the results con- 
cerning the mixing of operators are valid in general. Note also that our considerations 
apply equally well to moments of distribution amplitudes. 

Let us fix the notations used in our perturbative calculations. We work in Euclidean 
space and use the Wilson gauge action and Wilson fermions such that the total action is 
given by 

olatt q i C 

>->W — "->W,F + "JW,G • 
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The fermionic part <Sw,f for dimensionful massless fermion fields ip(x) has the form 
<Sw,f = Aa?r s y^ j ip(x)ip(x) 

X 

a 3 



x,ii 



where a is the lattice spacing and the sums run over all lattice sites x and directions fi 
on the lattice (all other indices are suppressed). The link matrices U X)li are related to the 
gauge field Ap(x) by 

U X)fl = exp [igaA^x)] , A^x) = T c A c ^x) , 

where g is the bare gauge coupling and the T c are the generators of the SU(3) algebra. 
The Wilson parameter r can be chosen from the interval (0, 1]. The gauge action for the 
gluon field A^(x) is 



5w,g — 

9 



1 - -Tr {U X4LV + Ul^) 



with 



X,f±<V L 



U x ,iw U x ^U x+a ^ tU U x+a p ^Ul u . 



In the perturbative calculation the investigated operators are sandwiched between off- 
shell quark states. We shall denote the momentum of the incoming quark by p and that 
of the outgoing quark by p' . Our calculations are performed in Feynman gauge, the final 
numbers will be presented for the Wilson parameter r = 1. 



2 Operators and mixing 



2.1 Renormalisation and mixing in the one-loop approximation 

In this section we discuss renormalisation and mixing in general terms, on the lattice as 
well as in the continuum. 

Let T^p', p, /i, gn, e) (j = 1,2, ...,N) be the dimensionally regularised amputated 
vertex functions of N mixing operators Oj of the same dimension calculated in 4 — 2e 
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dimensions. The corresponding Born terms are denoted by Tj n (p',p). The operators 
potentially contributing to the mixing have to satisfy certain symmetry requirements. 
They should transform identically according to a given irreducible representation of 0(4) 
or H(4), respectively, and they should have the same charge conjugation parity. 

The (dimensionless) renormalised coupling constant g R is related to the (dimensionful) 
bare coupling constant g by 

£ = ATV (1 + GOT 1 )) , (6) 

where \i is the renormalisation scale. In one-loop perturbation theory we get results of 
the form 



T?(p',p,^g R ,e) = T^(p>,p) 



+ 9r 



f>$ (j - + H**) - in ^±jPj rt or V,p) + W,V) 



0{g R ) 



(7) 

where 7^ = 0.5772 ... is Euler's constant. As usual, contributions which vanish for e — > 
have been omitted. In what follows, we systematically omit all contributions 0(g R ). In 
the MS scheme the renormalised vertex functions are then given by 

Tf(p',p,^g R )=Tf°™(p',p) 



9l 



N 



k=l 



Afi 2 



In the absence of mixing with lower-dimensional operators the vertex functions regu- 
larised on a lattice can be written as 



Tf(p',p,a,g R ) = rf ™(p>,p) 



+ 9r 



N 



£77 fe "(-I)" In 



ar(p' + p) x 



rf rn (p',p) + /;(p',p) 



k=l 



(9) 



up to terms vanishing as a — > 0. There should be an iV x iV matrix ( such that the relation 
between the bare lattice vertex functions and the MS renormalised vertex functions can 
be written as 



N 



rf (p',P, H, 9r) = ( 6 i k + g R^ k ) T k(p'iPi a > 9r) ■ 



(10) 



k=l 



So we should have 

^R(„> ~ g 
N 



Tf(p' 1 p J ^g R ) = Tf°™(p' 1 p) 

a 2 (p' +p) 2 



+ 9r 



k=l ^ 
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r B k mn (p', P ) + fHp',p) 



'11] 



Comparing with JSJ) we arrive at 



N 



£ fe* - t]> (aV 2 )) rn (p',p) + //(j/.p) - = . (12) 

k=l 

As this equation must hold for arbitrary momenta p', p, there should be constants cK 
such that 

ffip\p) - w,p) = £^rrv,p) • ( 13 ) 

fc=i 

This fixes the matrix Q^: 

<Jk = lpn(aV)-cJ k . (14) 

If the coefficients cj k are uniquely determined, exactly all N operators mix to one-loop 
accuracy. It might happen that certain operators contribute to mixing only in higher 
orders. In that case N is decreased. If Eq. (JT3J) cannot be satisfied, this shows that at 
least one mixing operator has been overlooked. 

Mixing with lower- dimensional operators leads to the appearance of additional terms 
on the r.h.s. of Eq. (jHJ). For example, consider the case of a single operator which mixes 
with Oj and whose dimension is one unit smaller. (Such a case will appear in our appli- 
cations.) Then we get instead of © 

rt(p',P,a,^) = rf m (p',p) 



+gi 



(15) 



itll ■ (-1) • m QV 4 +P) V r V,P) + V B -V,P) + f?(p', P ) 

(Jj 

,k=l 



where T BoTn (p' ,p) is the Born term of the additional lower-dimensional operator. The 
1/a contribution has to be subtracted from Tj(p',p, a, gn) before the connection with 
Tf{p',p,jjL,g^) can be established, i.e. in (JTUJ) T\{p',p, a, g R ) has to be replaced by 

T L k (p',p,a,g R )-^c k T^(p',p). (16) 
a 

Then Eq. (II 3|) is obtained as before. 

In order to compute the matrix Zjk of renormalisation and mixing coefficients we note 
that the connection between the bare lattice vertex functions and the MS renormalised 
vertex functions can be written as 

N 

Tf(p',p,fi,g R ) = Z^J2 Z ^(p',p,a,g R ) (17) 

fc=i 
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with the quark wave function renormalisation constant 2 (In the presence of mix- 

ing with a lower- dimensional operator, r%(p',p,a,gn) has again to be replaced by the 
subtracted expression |T6|) ). Comparison with (fTUj) yields 

Z^Zjk = $jk + OkQu ■ (18) 



The quark wave function renormalisation constant is calculated from the quark prop- 
agator. Write the lattice regularised inverse quark propagator in Feynman gauge (after 
subtracting a linearly diverging contribution ~ 1/ct) as 



S L \p) = i i> 



(19) 



where &l — 11.8524 (see e.g. ^H|) for Wilson fermions. In dimensional regular isat ion the 
inverse quark propagator is given by 



skc F 



167T 2 



S- D \p) = i 

and we get for the MS renormalised propagator 

S-\p) = if 

Defining such that 
we have 



1 p 2 
--JE + hi(47r) - In ^ 



(20) 



I 9 2 rC f J 




16vr 2 < 





S^(p) = z^ l s- L \p) 

^ = l-#f (ln(aV) + l + a L ). 



(21) 

(22) 
(23) 



With the help of this result we get from ()14|) and (jl8j) for the matrix of the renormalisation 
and mixing coefficients 



167T 2 



In (a 2 /x 2 ) 



;jk 



^16^ (1 + ^ 



(24) 



The basic computational task is thus to find the functions fj(p',p) and ff{j>'-,p) in 
Eqs. (J7J) and (JHJ) (or rather their difference). Then we can compute the coefficients cj k 
from Eq. (fT3*j) and use Eq. (j24j) to get the desired renormalisation and mixing coefficients. 
If we are only interested in the renormalisation of one particular operator, corresponding 
to j — 1 say, it is sufficient to restrict the calculations to the case j = 1. 

2 Often the quark wave function renormalisation constant is defined as 1/Z^,. 
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2.2 Contributing operators for special representations 



Let us introduce the self-explaining notations for operators with covariant and external 
ordinary derivatives 



C£ = -\d„^D^) , (25) 

oil = -^a(^) 



1 -r 



and 



1 



O 5 '™ = —dvUwsDuiP) , (26) 



1 



4 



0%™ = -- A d v d^i, W if) 



4 

as well as the lower-dimensional operators 

Og* = ~\$h» 7,] 5. V , OjL = --d u 7,M • (27) 

For completeness we include here also operators contributing to non-forward transversity 
matrix elements with two derivatives which were not considered in Eqs. (0) and (0J): 

= " jfe, 7,] V , 0££ = (fe, 7*4» • (28) 

As short-hand notations we use in the following (cf. |12j ) 

(C^l U2U3 + Cl/i 1/3^2 4~ U2UlU3 + C^ 2 ^3^1 ^^3^1 ^2 Ou 3 U 2 Ui) j (29) 
^1^2^3 OuiU3U2 ~l~ ^1/3 1/2 2 U2USUl + 2 C^j,-^ , (30) 

V\V$VZ ^^1^3^2 C^3^iI/2 ^31/21/1 • (31) 

First we consider the operator 

Cfi?4} = -^7{iB 2 B 4} ^. (32) 





1 


^{i/l 1/2^3} 


6 


Cm ii 

^|| ^1^2 ^3 || 


= o t 


0({^1^2^3>> 


= o, 



Its charge conjugation parity is C = — 1 and it is a member of an irreducible multiplet 
of operators transforming according to the representation of H(A) [T2|. Here rjp 
denotes an irreducible representation of H(4) with dimension I, and k = 1,2, . . . labels 
inequivalent representations of the same dimension. This operator can only mix with 

0{?24 } = -\d {2 d 4 (fayl,) . (33) 

Next we examine the operator 

Ox = OfA - \ (Of£ } + 0™ } ) . (34) 

It has already been used in lattice computations of forward hadronic matrix elements, 
because in this case it suffers only from rather mild mixing problems. It belongs to the 
representation t± with C — — 1. 

Taking into account also external ordinary derivatives, one finds the following opera- 
tors which transform identically and could therefore mix with 

2 = C{fl4} — ^ (^{224} + ^{334}) > 

n - n DD 1 (n DD + n DD \ 

U 3, - U ((114)) ~ 2 l U «224)> + U «334»J 
1 

'«114» - ^ ^«224» T- ^«334», 



3. 



04 — Cf{ll4)) - 7^ (^f{224)> + Cf(334») > (35) 
^5 — * — 'l (23.3| | ' 

U 6 - U ({213)) 
5,DD 
l|213|| 



o 7 = o 



and the lower- dimensional operator 



Os = O d 411 - \ {0% 2 + 0? 33 ) . (36) 



As an "axial" analogue of (J32|) we consider the operator 



5 $° = -^ {1 D 2 D 4y ^, (37) 



3 The charge conjugation parities C of and O5 coincide since O7 is antisymmetric in the indices of 
the covariant derivatives. 
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which can mix with 

°{i2% = ~\d { 2d 4 (^7i}75^) • (38) 
(j3ZD and (jSHD belong to r 3 (4) with C = +1. 

Similarly, we have as a counterpart of (|34|) the operator 

°! = O&S? - 5 (OJS? + OS?) ■ (39) 

Its charge conjugation parity is C = +1 and it is a member of an irreducible multiplet of 
operators transforming according the representation t% of H(4) [T2*] . 

The following operators with identical transformation behaviour could potentially mix 
with flUD: 



5 



o 



o\ = o 



,5,9a 

{114} 


5< 


s^i5,dd , x/-\5,aa A 

^{224} °{334} ^ ' 


,5,DD 
«114» " 


1 

" 2 


\ U ((224)) + 


^«334» 


,5,99 
«H4» " 


1 

" 2 


/ ^,5,99 
^((224» ■+■ 


sr\5,dd 
^«334» 


3D 








1 1 213 1 1 > 








3D 








((213)} > 








,DD 








||213|| 









^,5 _ /n5,9a _ ^ / ^,5,99 ^,5,99 \ /^x 
^4 — ^7(114^ o I ^((224)) ^ ^((334^ i ' 

e> 5 5 = o 
ojj = o 
0^ = 

and the lower- dimensional operator 

Ol = Og, - 20? 31 - Of 32 . (41) 

For the operators (|28|) we consider the representations r^ 3 ' 1 , Tg 3 " 1 and with C = 
—1 [T3] . In the case of r 2 we choose the representative operator 

Of = , (42) 

which may mix with 



Ol = • (43) 



For we take 



QT _ _ Q T,DD Q T,DD _ T,DD Q T,DD _ 2 qT,DD (u) 
W 3 — ^1{133} ^ ^1{144} ^2{233} ^ ^2{244} ZL/ 3{344} V^J 

11 



mixing with 

'l{133} ^ ^1{144} ^2{233} ^2(244} iv ^3{344} ' 



n T - -n T ' dd -u n T ' 9d - n T >" + n T > dd - on T ' 99 (a^\ 

^4 — ^111331 ^111441- ^21233) ^2|244> ZL/ 3l344l ■ 



The operator 

'l3{32} "T" *^23{31} ^14{42} ^24{41} 



n T — n T ' DD _l n T ' DD _ n T > DD _ n T ' DD (ar\ 

U 5 — U 13V32\ ^231311 L/ 14|42j ^24141) l 40 J 



belonging to the representation Tg 6 '* also mixes with only one additional operator: 

^6 — LV 13{32} LV 23{31} Ly i4{42} Ly 24{41} • V 11 ) 

Other representations show a more complicated mixing behaviour and will be not consid- 
ered here. 



3 One-loop calculation of non-forward matrix ele- 
ments 

3.1 Calculational technique 

We calculate the matrix elements of the operators in one-loop lattice perturbation theory 
in the infinite volume limit following Kawai et al. The computation is performed 

symbolically adopting and significantly extending a Mathematica program package devel- 
oped originally for the case of moments of structure functions using Wilson ^3] , clover 
and overlap fermions [T7j . 

Using that approach we have full analytic control over pole cancellation. The Lorentz 
index structure of the matrix elements is left completely free, so that we are able to 
construct all representations of the hypercubic group for the second moments in the non- 
forward case. The analytic derivation is disconnected from the numerical calculation of 
finite lattice integrals. These integrals are calculated to a high accuracy and are given as 
look-up tables. 

To be more precise, we recapitulate in short the strategy Ej for calculating the 
diagrams adapted to our case of two distinct external momenta p 7^ p'. Basically we have 
to evaluate a typical lattice integral of the form 

r 7T / a (fik 

Zn-HnMiP) = / 7^u^w-**»( a >P'>P> fc ) ( 48 ) 

J-Tv/a (A) 
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where the integration is performed over the first Brillouin zone. The numerator of the 
kernel K. fll ... IMn (a,p',p,k) is a polynomial in sines and cosines of the lattice momenta k, 
the denominator contains the denominators of lattice quark and gluon propagators. Such 
an integral X is calculated by splitting it into two parts 

1 = 1 +(1-1). (49) 

Here 1 denotes the Taylor expansion of the original integral 1 in the external momenta 

l(a,p',p) = X(a,0,0) 

/ / dl(a,p',p) dl(a,p',p) \ , . 

+ 2^\p« dg „_ n + P«— £T — + ■■■ ( 5 °) 



P '=p: 



=o dp a 



where the order of the expansion is given by the degree of ultraviolet (UV) divergence of 1. 
As a consequence, the difference 1—1 is UV finite and can be calculated in the (Euclidean) 
continuum, taking the limit a — > 0. The original UV poles appear now as infrared (IR) 
poles in the Taylor expansion and are regularised using dimensional regularisation with 
d > 4. Since in that case l\ a ->o vanishes, we have (1 — X)| a _ )0 = X| a ^ = l cont (p',p) and 
the second contribution to the final lattice result is just a one-loop continuum calculation 
in dimensional regularisation. 

The first (Taylor expanded) part is calculated in d = 4 — 2e dimensions at finite a, 
the poles in e analytically cancel those of the second part X cont . Note that the first part 
depends on the external momenta only via the Taylor expansion, therefore the lattice 
integrals (the expansion coefficients) are just numbers (independent of p and p'). 

In the Euclidean continuum part l cont (p',p) the finite contribution depends on the 
(different) external momenta in a complicated manner. As in Minkowski space, the finite 
part of three-point functions with three different propagators is difficult to represent in a 
simple compact analytic form (see, e.g., [IBj)- in parametrising that contribution we have 
chosen the following form 

^° n V>>) = - € ^i-^(p',p) + B^...^{p',p,i,j,k,m)T(i,j,k,m,p',p) . (51) 

i ,j,k,m 

In (pTTj) the functions A^...^ (p', p) (known analytically) and B^...^ (p', p, i, j, k, m) contain 
the general index structure assigned to the external momenta p and p' . The symbols T 
stand for the remaining integrals (finite for p ^ p') over the Feynman parameters 4 : 

F(i,j,k,m,p',p)= f dx f X dy x* j? (Q 2 (x, y,p',p)) k In™ Q * (x ' v ; P '' p) (52) 



One Feynman parameter integration can be performed analytically, leading, however, to lengthy 
expressions for the individual integrals. We have checked that for our continuum results we recover the 
known results of the forward case. 
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with 

Q 2 (x,y,p',p) = p 2 x(l -x) + p' 2 y(l - y) -2p-p'xy. 

Only few combinations of the integers i,j,k,m are actually needed in the "continuum 
part" of the one-loop calculation. But compact expressions for the !F(i,j,k,m,p',p) do 
not seem to exist. 

As we have shown in Section ITT] in order to calculate the finite renormalisation matrix 
converting the lattice result to the MS scheme we have to find the difference between the 
lattice and the continuum one- loop contribution ff{jp',p) — fj{p',p), and similarly for the 
wave function renormalisation. Therefore, we do not need the explicit form of the finite 
continuum one- loop contribution (which is essentially given by the .F's), since it cancels 
exactly the part coming from the second term of (J49|) . As a consequence, we do not 
present jF-integrals here. However, choosing other schemes, the finite contributions to 
^Mi - Mn(p''P) cont m ight be needed. 

3.2 One-loop renormalisation matrices for the chosen represen- 
tations 

To calculate the one-loop contributions for the above operators between quark states 
we have to take into account self energy and amputated Green function diagrams. The 
quark self energy diagrams (contributing to the quark wave function renormalisation) are 
shown in Fig. [T] Denoting the operator insertions by a black dot, the vertex and tadpole 



k 




p p + k p p 



Figure 1: One- loop diagrams contributing to the quark self-energy. 

diagrams (Fig. EJ) as well as the "cockscomb" diagrams (Fig. |3J) contribute in one-loop 
order to the amputated Green functions. 

We present the results for the Z factors of the operators discussed in Section 12.21 in 
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Figure 2: One-loop vertex and tadpole diagrams. 




Figure 3: One-loop cockscomb diagrams. 

the generic form 

z * = - w ln(aV) + C S } ) (53) 

where (cf. (El) 



167T 2 v ( m ) 167T 2 , 

7jfc = - -^-7j* , c }fc = M 1 + + ~c^ c jk ■ ( 54 ) 



The superscript (m) with m — 1,11 distinguishes the realisations I and II of the covariant 
derivatives, which are explained in the Appendix. 
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3.2.1 Of° i} (r?\c = -l) 

For this representation we find mixing between 0^f 4} ® and Off 24} The corre- 

sponding 2 x 2-mixing matrices are 

/ 25 _5 \ 

7i* =(' ]). (55) 

_ / -11.5632 0.0241 \ 
C ^ fc ~ V 20.6178 J ■ W 

The anomalous dimension matrix was calculated earlier [TTJ. The element en = —11.5632 
in the matrix (|56|) is known from our previous calculation in the forward case ^3] and C22 
comes from the renormalisation of the local vector current ^H]- The matrix c^.' 7 ^ shows 



a very small mixing between the operators OP^y and 09® 24 y. Thus it may be justified to 
neglect the mixing in practical applications. 



3.2.2 C 1 (r 1 (8) ,C = -l) 



We consider the mixing of operators having the same dimension first. The relevant oper- 
ators have been defined in Section 12.21 To one-loop accuracy the operator O-j in Eq. (j35j) 
does not contribute and we have to consider the following mixing set: 

{Ol, 02, 03, 04, 05, 06}. 



The anomalous dimension matrix is 



\ 



25 
6 









5 
6 












7 
6 










5 
6 










1 


2 

2 






3 
2 


-2 

2 
3 



and the finite part of the mixing matrix is given by 

/ -12.1274 -2.7367 0.3685 0.9934 

20.6178 

3.3060 18.1841 -14.8516 -4.3023 

20.6178 

3.2644 

\ 3.2644 



-jk 



0.0156 


-0.9285 


0.3501 
0.0050 



0.1498 \ 


0.7380 


0.0149 
0.3600 / 



(57) 



(58 
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or 



c jk 



( 


-12.1274 


1.4913 


0.3685 


-0.4160 


0.0156 


0.1498 \ 







20.6178 
















3.3060 


-8.01456 


-14.8516 


4.3023 


-0.9285 


0.7380 













20.6178 













3.2644 








0.3501 


0.0149 


\ 





3.2644 








0.0050 


0.3600 / 



(59) 



The matrices c^.' show a sizeable mixing of the operator 0\ with other operators, 
especially O2 containing two external ordinary derivatives. This mixing becomes irrelevant 
in the forward case, where the matrix element of O2 vanishes. On the other hand, the 
mixing between the operators 0\ and O3 is already known from our previous calculations 
in the forward case [TB] . 

There is also a possible mixing between 0\ and the lower- dimensional operator 8 (|36|). 
Indeed, we find in the one-loop approximation that the vertex function of 0\ contains a 
term oc 1/a: 



01 



-part 



(-0.5177) ^0 8 Born 



(60) 



This mixing leads to a contribution which diverges like the inverse lattice spacing in the 
continuum limit. The operator Os has to be subtracted non-perturbatively from the 
operator 0\ which might be a difficult task in simulations. 



3.2.3 Olf 2 °(rt\C = +l) 

For this representation we find mixing between and C^f} ■ The anomalous dimen- 

sion matrix is given by (|55|). the finite contributions are collected in the matrix 

(/,//) _ / -12.1171 0.1667 
C * ~ V 15.7963 I ■ 

As in the case of the operator O?^, above, the mixing is rather small, and the diagonal 
matrix elements en and C22 can be inferred from previous work jT 
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3.2.4 0\(rf\C = +1) 



First we discuss the mixing of operators of the same dimension. The set of contributing 
operators is found to be 

{olololololot}. 

As in the case of the operator Oi, one operator - here Oj - does not contribute to mixing 
in one-loop order. The finite contributions for the two lattice derivatives are 



and 



( -12.8609 


3.4220 



V o 

/ -12.8609 


3.4220 







-2.0653 
15.7963 
15.8207 


-8.9124 
-8.9124 



1.4894 
15.7963 
-7.3002 


-8.9124 
-8.9124 



0.3490 


-15.3592 






0.3490 


15.3592 






0.8538 


-5.1639 
15.7963 





-0.3311 


2.5431 
15.7963 





0.0511 


0.1701 


0.9597 
-0.3199 



0.0511 


0.1701 


0.9597 
-0.3199 



0.0594 \ 


-0.9431 


-0.9597 
0.3199 / 

0.0594 \ 


-0.9431 


-0.9597 
0.3199 / 



(62) 



(63) 



The anomalous dimension matrix is the same as for the operators without 75, see (JSTj) . 
Again, some of the mixing coefficients may be non- negligible. The mixing between 0\ 



and 0\ is visible even in the forward case. 



We also find mixing with a lower- dimensional operator, in this case it is the operator 
0\ (|41j) . The corresponding contribution in the vertex function of 0\ reads 



-part 



(64) 



3.2.5 0\ (r 2 (3) , C = -1) 

For the mixing between the operators 0\ f!42j) and 0\ (}4*3*|) the matrix of anomalous 
dimensions is given by 

7i* = ( I ~f ) , (65) 
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and the finite contributions are 



n n) _ ( -11.5483 0.2189 



3.2.6 Ol{rf\C = -l) 



The mixing between the operators Oj f|44|) and Oj (|45j) leads to the same anomalous 
dimension matrix (|65|) as the previous case. For the finite pieces we obtain 

(/,//) _ / -11.8688 0.2753 \ 
C i k ~ I 17.0181 r [b<) 



3.2.7 OUt?\C = -1) 

For the operators 0\ f|46|) and Oj ()47|) we find the finite mixing contributions 

(i,n) _ ( -11.7477 0.2380 \ 

* \ 17.0181 ) 1 

with the identical anomalous dimension matrix (|63j) . 



4 Summary 



Within the framework of lattice QCD with Wilson fermions and Wilson's plaquette action 
for the gauge fields we have calculated the one-loop quark matrix elements of operators 
needed for the second moments of GPDs and meson distribution amplitudes. From these 
we have determined the matrices of renormalisation and mixing coefficients in the MS- 
scheme. The operators contributing to the considered moments in the continuum contain 
two derivatives and are hence of dimension 5. On the lattice we should take operators with 
definite C-parity belonging to multiplets which transform irreducibly under the hypercu- 
bic group. We have selected some operators commonly used in numerical simulations of 
ordinary parton distributions. While these distributions are extracted from forward ma- 
trix elements, we have to consider non-forward matrix elements in order to study GPDs. 
This circumstance enlarges the set of potentially mixing operators, because also operators 
with external ordinary derivatives have to be taken into account. 
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For the representations (r 2 ,C = — 1) and (r 3 , C = +1) there are only two multiplets 
of operators of dimension < 5. Actually, they are all of dimension 5, and the mixing 
pattern is the same as in the continuum. The mixing coefficients turn out to be quite 
small. 

In the case of the representations (rf , C = —1) and (r 2 ,C = +1) we have eight 
multiplets of operators of dimension < 5. In the one- loop approximation one of them 
cannot yet contribute, but the remaining seven multiplets do actually show up. Six of 
them are of dimension 5, but there is also one multiplet of operators of dimension 4. 
For these representations some of the mixing coefficients are larger than in the above 
cases. In particular the mixing with the lower-dimensional operators, which will probably 
receive sizeable non-perturbative corrections, might lead to difficulties in the numerical 
simulations: 1/a effects are hard to get under control. 

In the operators just mentioned the Dirac matrix is either 7^ or 7^75. In connection 
with (generalised) transversity distributions operators with [7^, 7„] are of interest. Here we 
restricted ourselves to cases where only mixing with the external derivative counterparts 
is possible and found rather small mixing coefficients. 

The results presented in this paper were obtained for Wilson fermions. However, nu- 
merical simulations will be performed with clover or overlap fermions. Therefore, the 
numbers given in the preceding section can only serve as first hints at the problems which 
will occur in these calculations, especially the mixing problem. The sets of mixing oper- 
ators, on the other hand, follow from symmetry arguments and are therefore universally 
valid, although additional symmetries may lead to further restrictions, e.g. in the case of 
overlap fermions. It should also be noted that the overall normalisations of the opera- 
tors are somewhat arbitrary. So some care has to be exercised when using our results in 
concrete applications. 
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Appendix: Feynman rules 



We need Feynman rules for the lattice operators (Jzojl . (126(1 and (|28|) . Denoting the mo- 
mentum of the outgoing quark by p', that of the incoming quark by p and the incoming 
gluon momenta by ki, we use the following Fourier decomposition for the quark and gauge 
fields, 



d A p 

d 4 p' 
d 4 kj 

(27T) 4 



i>(p')e- ip '- x 



(69) 



A^ki) e 



iki-(x+afl/2) 



where the momenta are restricted to the first Brillouin zone. 

We use the standard realisation of the covariant derivatives acting to the right and to 
the left 



Dp 4>{x) 



2a L 
1 

2a L 



U Xt(t ip(x + afi) - Ul_ a ^ ij)(x - aft) 



ip(x + afi) Ul — ^jj(x — aft) U x . 



ap : p 



The external ordinary derivative is taken as 

dp {4> - ■ ■ V>) (x) = - [($ ■ ■ ■ ip) (x + afi) - • • • if) (x)] 



(70) 



(71) 



There are two convenient ways to define the corresponding operators in momentum 
space with non-zero momentum transfer q. One possibility is to "act" with q at the lattice 
point x. For the case of one covariant derivative this leads to (setting the Dirac matrix 
in the operator equal to the unit matrix for simplicity) 



Dp ^} {I \q) = Y,{4> Dp 4>)(x) e l 

X 

= Ya S [^ X ) U ^( X + a ^ - fa + *fiUijK*j\ ¥ q ' X + e*^] . (72) 

X 

Alternatively, the momentum transfer can be applied to the "position center" (for the 
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case of one covariant derivative D) as follows 



x—an,u 

which leads to 

(H) 



*i> ij) (q) 

X 

This realisation might be more suitable for numerical simulations. Note that 

D[i tp) (q) = cos DfM tp) (q) . 
For the external ordinary derivative we use 



a 



= ~~ sin ~Y (V"-"0)(9)- 

The corresponding Feynman rules up to 0(g 2 ) derived from expressions based 
and (fTHJ) are marked by the superscripts (I) and (II), respectively. 

O(g ) 



n (H),DD 




1 c 
— sin - 

a 1 


n (I),DD 

flUUI 


a(p - p') 

= cos - 

2 


v a 
- cos — 


n (II),8D 


= *P(p'hti*P(p) 


1 I 

— sm - 

a z 


n (I),dD 

^ flUUJ 


a(p- p') 

= cos 

2 


« wj/) 


Qdd 




1 c 
_ sm- 
a 2 



sm 



£_ /n(II),DD 
2 jj,VU) 1 

1 . a{p-p') v . a(p + p') u 
sm 



1 a{p-p') v . a(p-p'), 
sm 



22 



0(9) 



(II),DD 



n {I),DD 
n {II),dD 

^ (1VU) 

n (i),dD 



g^2vp(p')>y^A a (k 1 )'ifj(p) cos 

C 

1 



a(p + p') a 



. a{p + v'-k\)u , j- a(p + p' + k 1 ) l 
Ova sin V d WCT sin 



COS 



a{p—p' + ki) v a(p — p' + kx) 



cos 



w q(II),DD 



(77) 



gfp{ph^ u ,(ki)i/j(p)-sm cos 



cos 



a(p - p' + fci)^ ^(//j.ar, 



0(g 2 ) for the tadpole case k 2 = —k\ 

Here we restrict ourselves to the tadpole case k 2 = —k\ (cf. Fig.[2J) as this is all we need 
in this paper. For general gluon momenta the Feynman rules are much more complicated. 

<T,T 

ov c a(p + p' + k x ) v a(p + p' + ki) u 
x < Zo TV o mil cos cos 



fx x i x x \ ' a (P + P')» ■ a (P + P% 

- {OruOau + TW d au} ) Sin Sin 



0^ D = cos cos ^^^^ 

= -^(^7,^(^)^(-^)^(p)sin^^sin a(P + ^ 



2 7 VJ ' ,p v ' v ' T " ' 2 2 

0Wf = cos a(p ~ P% 0^ 8D . 

flVU) ry (IV UJ 

The Feynman rules for the operators with different Dirac matrices are obtained by the 
obvious replacements. Note that the forward case is realised for 0(g°) and 0(g 2 ) (in the 
tadpole case) by p = p', for 0(g) we have to take p' — p + k±. 
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